# Simulate Fall 3NO female plaice population
## See "imitate_plaice_fall_3NO_female.R" file for details on parameter choices
library(SimSurvey)
set.seed(889)
pop <- sim_abundance(ages = 1:26,
years = 1:20,
R = sim_R(log_mean = log(100000000),
log_sd = 0.7,
random_walk = FALSE),
Z = sim_Z(log_mean = log(0.15),
log_sd = 0.5,
phi_age = 0.9,
phi_year = 0.5),
N0 = sim_N0(N0 = "exp", plot = FALSE),
growth = sim_vonB(Linf = 69.63, L0 = 3, # Fitted for female growth
K = 0.09, log_sd = 0.1,
length_group = 2, digits = 0)) %>%
sim_distribution(grid = make_grid(x_range = c(-184, 184),
y_range = c(-184, 184),
res = c(3.5, 3.5),
shelf_depth = 60,
shelf_width = 170,
depth_range = c(0, 1600),
n_div = 2,
strat_breaks = seq(0, 1600, by = 65),
strat_splits = 4,
method = "bezier"),
ays_covar = sim_ays_covar(sd = 1,
range = 800,
phi_age = 0.9,
phi_year = 0.9,
group_ages = 20:26),
depth_par = sim_parabola(mu = log(75),
sigma = 0.1,
sigma_right = 0.55, log_space = TRUE))
## Test a series of surveys
## Simulate surveys and compare stratified estimates to the true index
## Include baseline, increase and decrease of 80% and 50%; except for age, only 50%
surveys <- expand_surveys(set_den = c(0.5, 0.8, 1, 1.2, 1.5) / 1000,
lengths_cap = c(75, 120, 150, 180, 225),
ages_cap = c(10, 20, 30))
surveys[surveys$set_den == 0.001 &
surveys$lengths_cap == 150 &
surveys$ages_cap == 20, ] ## survey 38 ~current protocol
sim <- test_surveys(pop,
surveys = surveys,
keep_details = 38,
n_sims = 5,
n_loops = 200,
cores = 7,
q = sim_logistic(k = 2, x0 = 2),
export_dir = "C:/Users/FITZSIMMONSM/Documents/SimSurvey/SimSurvey/analysis/plaice_female_sim_exports/2021-01-11_age_clust_test")
# sim <- resume_test(export_dir = "analysis/cod_sim_exports/2018-10-26_age_clust_test")
## Visualize Results
load("C:/Users/FITZSIMMONSM/Documents/SimSurvey/SimSurvey/analysis/plaice_female_sim_exports/2021-01-11_age_clust_test/test_output.RData")
plot_total_strat_fan(sim)
plot_length_strat_fan(sim, surveys = 1:75, years = 1:20, lengths = 0.5:144.5)
plot_age_strat_fan(sim, surveys = 1:75, years = 10, ages = 1:26)
plot_error_surface(sim)
plot_survey_rank(sim)
## Test alternate survey with strat specific age sampling and age-length-keys
surveys <- expand_surveys(set_den = c(1, 1.2) / 1000,
lengths_cap = c(75, 120, 150, 180, 225),
ages_cap = c(10, 20, 30))
sim <- test_surveys(pop,
surveys = surveys,
keep_details = 1,
n_sims = 5,
n_loops = 200,
cores = 7,
q = sim_logistic(k = 2, x0 = 2),
export_dir = "C:/Users/FITZSIMMONSM/Documents/SimSurvey/SimSurvey/analysis/plaice_female_sim_exports/2021-01-14_set_alk",
age_length_group = 2,
age_space_group = "set",
alk_scale = "set")
# sim <- resume_test(export_dir = "analysis/cod_sim_exports/2020-02-10_set_alk")
## TODO: perhaps run sim_survey_parallel and get fan plots working for that ouptut?
## Visualize Results
load("C:/Users/FITZSIMMONSM/Documents/SimSurvey/SimSurvey/analysis/plaice_female_sim_exports/2021-01-14_set_alk/test_output.RData")
plot_age_strat_fan(sim, surveys = 2, ages = 1:10, years = 3)
plot_survey_rank(sim, which_strat = "length")
## Simulate same distribution across ages --------------------------------------
## Same as above cod-like simulation, except spatial distribution is the same
## across all ages
rm(pop)
rm(sim)
gc()
set.seed(889)
pop <- sim_abundance(ages = 1:26,
years = 1:20,
R = sim_R(log_mean = log(100000000),
log_sd = 0.7,
random_walk = FALSE),
Z = sim_Z(log_mean = log(0.15),
log_sd = 0.5,
phi_age = 0.9,
phi_year = 0.5),
N0 = sim_N0(N0 = "exp", plot = FALSE),
growth = sim_vonB(Linf = 69.63, L0 = 3, # Fitted for female growth
K = 0.09, log_sd = 0.1,
length_group = 2, digits = 0)) %>%
sim_distribution(grid,
ays_covar = sim_ays_covar(sd = 1,
range = 800,
phi_age = 0.9,
phi_year = 0.9,
group_ages = 20:26),
depth_par = sim_parabola(mu = log(75),
sigma = 0.1,
sigma_right = 0.55, log_space = TRUE))
## Test a series of surveys
## Simulate surveys and compare stratified estimates to the true index
surveys <- expand_surveys(set_den = c(0.5, 1, 2, 5, 10) / 1000,
lengths_cap = c(5, 10, 20, 50, 150, 500, 1000),
ages_cap = c(2, 5, 10, 20, 50))
surveys[surveys$set_den == 0.001 &
surveys$lengths_cap == 150 &
surveys$ages_cap == 20, ] ## survey 127 ~ roughly current protocol
sim <- test_surveys(pop,
surveys = surveys,
keep_details = 127,
n_sims = 5,
n_loops = 200,
cores = 7,
q = sim_logistic(k = 2, x0 = 2),
export = "analysis/cod_sim_exports/2018-10-28_no_age_clust_test")
# sim <- resume_test(export_dir = "analysis/cod_sim_exports/2018-10-28_no_age_clust_test")
# ## visualize results
# load("analysis/cod_sim_exports/2018-10-28_no_age_clust_test/test_output.RData")
# vis_sim(sim)
## Test survey with set-specific age-length-keys --------------------------------------
library(SimSurvey)
set.seed(889)
pop <- sim_abundance() %>%
sim_distribution()
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.